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Abstract 

We present a search for the decays — > ^"'"f^ and — > e'^i'e in a 253 fb~^ 
data sample collected at the T(4S') resonance with the Belle detector at the KEKB 
asymmetric-energy B factory. We find no significant evidence for a signal and set 
90% confidence level upper limits of B{B+ ^i+u^) < 1.7 x 10"^ and B{B+ 
e+Ve) < 9.8 X 10-7. 
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1 Introduction 

The purely leptonic decay 5+ — > i^ug (charge conjugate states are implied 
throughout the paper) is of particular interest since it provides a direct mea- 
surement of the product of the magnitude of the Cabibbo-Kobayashi-Maskawa 
matrix element, \Vub\ [1], and the B meson decay constant, /b- In the Standard 
Model (SM) the branching fraction of the decay —>■ I'^vg is given as 

B(B* ^ I*.,) = { 1 - 4)'/I|V..I%. (1) 

where is the Fermi coupling constant, mi and are the charged lepton 
and B meson masses, and tb is the B^ lifetime. The expected branching 
fractions are (4.7 ± 0.7) x 10"^ for B+ fi+u^, and (1.1 ± 0.2) x 10"^^ for 
B^ e'^h'e assuming \Vub\ = (4.39 ± 0.33) x 10~^ determined from inclusive 
charmless semileptonic B decay data [2], tb = 1.643 ± 0.010 ps [2], and = 
0.216 ± 0.022 GeV obtained from lattice QCD calculations [3]. However, non- 
SM physics could yield larger branching fractions [4]. 

The Belle experiment has recently found the first evidence for the purely 
leptonic decay B^ —> t^Uj. [5]. The other purely leptonic B decays, B~^ — > 
yU"*"//^ and B^ — s> e'^Ue, have not yet been observed. The most stringent current 
upper limits for these modes are B{B^ — >■ yU^f^) < 6.6 x 10~^ [6] and B{B^ 
e+i/e) < 1-5 X 10"^ [7]. Preliminary limits of b{b+ /i+z/^) < 2.0 x 10"^ [8] 
and B{B^ e^Ve) < 7.9 x 10"^ [9] are also available from the Belle and 
BaBar collaborations, respectively. In this paper, we present a search for the 
decays B^ — > yU+i/^ and B^ — >• e+fg. 



2 Data Set and Experiment 

The data used in this analysis were collected with the Belle detector at the 
KEKB asymmetric-energy e+e~ collider [10]. The sample has an integrated 
luminosity of 253 fb~^ accumulated at the T(4S') resonance, at a center-of- 
mass (CM) energy of 10.58 GeV (on- resonance) , and 28.1 fb~^ accumulated at 
a CM energy 60MeV below the T(4S') resonance (off-resonance). 
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The Belle detector is a large solid-angle spectrometer based on a 1.5 T super- 
conducting solenoid magnet. Charged particle tracking and momentum mea- 
surements are made with a three-layer, double-sided silicon vertex detector 
(SVD) and a central drift chamber (CDC). Identification of charged hadrons 
is provided by a combination of three measurements: specific ionization loss 
(dE/dx) in the CDC, photon yield in the aerogel threshold Cherenkov counters 
(ACC), and time-of-fiight information from a cylindrical array of 128 scintilla- 
tion counters (TOF). Photons are detected in an electromagnetic calorimeter 
(ECL) system made of an array of 8736 CsI(Tl) crystals surrounding the TOF 
system. The ECL is also used for electron identification. Muons are identified 
by a resistive plate chamber system (KLM) located within the solenoid's ex- 
ternal return yoke. 

Particle identification for and /i^ is important for this analysis. Electron 
identification is based on the ratio of the cluster energy in the ECL to the track 
momentum from the CDC and the SVD (E/p), the dE/dx in the CDC, the 
position and shower shape of the cluster in the ECL and the response from the 
ACC. Muon identification is based on the hit positions and the depth of pen- 
etration into the KLM. The efficiency of electron identification is over 90 % in 
the momentum range of this analysis while the misidentification rate is below 
0.5 %. The muon identification efficiency is approximately 85 % in the momen- 
tum range of this analysis with a misidentification rate of approximately 1 %. 
The Belle detector is described in detail elsewhere [11]. 

A detailed Monte Carlo (MC) simulation, which fully describes the detector 
geometry and response based on GEANT [12], is used to estimate the sig- 
nal detection efficiency and to study the background. The MC samples are 
produced with the EvtGen event generator [13]. 



3 Signal Selection 

We search for events in which there is one well-identified lepton. The signal 
candidate muon or electron is required to pass tight particle identification 
requirements; lepton candidates that can be associated to other tracks in the 
event to reconstruct mesons or photon conversions are explicitly vetoed. 
Signal lepton candidates are required to originate from the interaction point 
(IP) and to have their polar angle, 6*^, formed by the lepton momentum with 
the detector axis (opposite to the direction of the positron beam), in the range 
—0.5 < cos^£ < 0.8 for muons and —0.50 < cos 6'^ < 0.85 for electrons. 

As i^i^e is a two-body decay, the lepton has a fixed momentum in the 

signal B meson (B^^^) rest frame, with pf equal to approximately half of the 
B meson mass, pf ~ 7713/2. The lepton momentum in the CM frame, p}, is 
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related to pf by 

^^-p;fl-^cos^,_sJ (2) 



where p^sig is the momentum of the B^^^ in the CM frame and cos^^.^sig 
represents the cosine of the angle between the directions of the signal lepton 
and -B*^'^ in the CM frame. 

Since the neutrino is not detected in the B'^^^ decay, we obtain p^sig by inclu- 
sive reconstruction of the companion B meson (5™™?) recoiling against B^^^. 
For the B^°^'^ reconstruction, we use all detected photons and charged tracks, 
except for the signal lepton candidate. Kg 7t^it~ and 7 e~^e~ decays 
are fully reconstructed, in order to correctly account for the momentum of 
tracks originating from vertices displaced from the IP. The missing momen- 
tum in the laboratory frame that is calculated by using all photons, charged 
tracks, and the signal lepton, is assigned to the neutrino. The quantity |p^Big | 

is approximately given by E'^^^^ — m\ ~ 0.32 GeV/c using the beam energy 
in the CM frame (-Ebeam), while cos6'£_^sig is related to the angle between 
the direction of the signal lepton and the momentum of B^°'^'^{9i_b^o^p) by 
cos6'£_5sig = — cos6'£_Bcomp. Thus, equation 2 can be expressed in terms of two 
measurable quantities, and cos^^-^comp: 

pf ^ pI{1 + 0.06 cos 6'£_Bcomp). (3) 



Signal candidates are selected using the kinematic variables Mbc = y -E'boam ~ \'Pb 
and Ai? = E*^ — -Ebcam? where and E*^ are the momentum and energy of 
^comp^ all variables being evaluated in the CM frame. Events in the fit region 
that satisfy 5.10 GeV/c^ < Mbc < 5.29 GeV/c^ and -0.8 (-1.0) GeV < < 
0.4 GeV for the muon (electron) mode are kept for further analysis. A more 
restricted signal region is also defined by the same requirements on AE' and 
by Mbc > 5.26 GeV/c^. 

The dominant background arises from the continuum process e^e" ^ qq {q = 
u, d, s, c) and semileptonic B meson decays, mostly into charm final states 
[B Xcd-v) with contributions from rare h u modes [B Xud-v)- The 
missing momentum of continuum events is often due to undetected particles 
that are outside the detector acceptance. In order to reduce such backgrounds, 
we require the transverse component of the missing momentum to be greater 
than 1.75GeV/c, and the cosine of the polar angle to be less than 0.84 (0.82) 
for the muon (electron) mode, respectively. 

Figure 1 shows the measured pf distributions after reconstruction of the com- 
panion B meson. We require 2.6GeV/c < p^ < 2.84GeV/c and 2.6GeV/c < 
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p» [GeV/c] p= [GeV/c] 

Fig. 1. Lepton momentum distributions in the rest frame of the signal B meson for 
the muon (left) and electron (right) mode after reconstruction of the companion 
B meson. Points show the on-resonance data, and solid histograms show the ex- 
pected background due to rare B — > Xuii^ decays (hatched, from MC); other BB 
events, principally B — > Xcii^ decays (cross-hatched, also from MC); and continuum 
events (light shaded, taken from scaled off-resonance data). The dashed histograms 
represent the signal as predicted by the MC with arbitrary normalization. 

pf < 2.8GeV/c for the signal candidates. This requirement removes most of 
the B — > Xcii' background. 

In order to further suppress the continuum background we exploit the event 
shape difference between continuum events, which are jet-like, and BB events, 
which tend to have a more spherical topology. Modified Fox- Wolfram mo- 
ments [14,15] are combined into a Fisher discriminant (F). Figure 2 shows the 
Fisher discriminant distributions for events that passed the cos 6*^ selection. 
We require F > 0.3 for the muon mode and F > for the electron mode, 
which retain approximately 51 % (60 %) of the signal in the signal region and 
remove approximately 99 % (95 %) of the continuum background in the signal 
region for the muon (electron) mode. 



We determine all selection criteria by maximizing Ns / y/Ns'+l^, where Ns is 
the number of signal events expected in the signal region computed assuming 
the SM branching fractions and Nb is the number of expected background 
events in the signal region from MC. 

After all selection criteria have been applied, the signal selection efficiency in 
the signal region is estimated to be = 2.18 ± 0.06% for the muon mode 
and ee = 2.39 ± 0.06 % for the electron mode. The efficiencies in the fit region 
are 3.15 ± 0.07 % for the muon mode and 3.86 ± 0.08 % for the electron mode. 
Figure 3 shows the pf distributions after all other selections have been applied. 
The background that remains in the signal region consists of approximately 
76 % continuum and 24% B ^ X^iu according to the off-resonance data and 
MC studies. The MC study also indicates that the latter exhibits a peaking 
Mbc distribution with significantly larger width than the signal distribution. 
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Fig. 2. Fisher discriminant distributions for the muon (left) and electron (right) 
mode after cosOg requirements have been applied. Points show the on-resonance 
data, and solid histograms show the expected background due to rare B — > JCy^£u 
decays (hatched, from MC); other BB events, principally B — > Xcii' decays 
(cross-hatched, also from MC); and continuum events (light shaded, taken from 
scaled off-resonance data). The dashed histograms represent the signal as predicted 
by the MC with arbitrary normalization. 




p= [GeV/c] p= [GeV/c] 

Fig. 3. pf distributions for the signal candidates. Points show the on-resonance data, 
and solid histograms show the expected background due to rare B — > X^j^z/ decays 
(hatched, from MC); other BB events, principally B — > Xcii^ decays (cross-hatched, 
also from MC); and continuum events (light shaded, taken from scaled off-resonance 
data). Dashed histograms are MC B ^ iv signals that are obtained by multiplying 
the SM expectations by a factor of 10 for the muon mode and by 5 x 10^ for the 
electron mode. The arrows show the signal regions. 

4 Signal Extraction 



The signal yields are extracted from unbinned maximum likelihood fits to Mbc 
distributions in the fit region. The signal Mbc distribution is parameterized 
by a Crystal Ball function [16] modeled on the signal MC. The background is 
described by ARGUS functions [17], with shape parameters determined from 
the off-resonance data for the continuum background and from MC for the BB 
background. We do not include the peaking component from B — > X„£z/ in 
this fit since an examination of AE sideband data indicates that the peaking 
contribution is negligible. The expected number of background events in the 
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Muon Mode Electron Mode 

Signal Efficiency (fit region) 
Signal Efficiency (signal region) 
Observed in Signal region [events] 
Expected background [events] 
Signal yield [events] 
Significance 
SM Prediction [events] 
Table"! 

Selection efficiencies, expected numbers of events for signal and background and fit 
results. 

signal region is estimated by fitting the Mbc distribution outside the signal 
region to a background shape determined from the BB MC and off-resonance 
data. The expected number of background events is 7.4 ± 1.0 (13.4 ± 1.4) 
for the muon (electron) mode. Use of the combination of the BB MC samples 
and off-resonance data instead of the on-resonance data gives similar expected 
number of background events. The Mbc distributions are used as probability 
density functions (PDF) to compute an extended likelihood function defined 
as follows: 

g-(ns+ni,) N 

C{ns, Uh) = — Ylirisfsii) + ^fe/blO) (4) 



3.15 ±0.07' 
2.18 ±0.06' 
12 

7.4 ± 1.0 
4.1 ±3.1 

1.3 
2.8 ±0.2 



7o 3.86 ± 0.08 % 
^0 2.39 ± 0.06 % 
15 

13.4 ± 1.4 
-1.8 ±3.3 

(7.3 ± 1.4) X 10" 



where Ug and rib represent the numbers of signal and background events in 
the fit region to be determined in the fit, iV is the number of observed events, 
fs and fb are the signal and background PDFs, respectively. The negative log 
likelihood function is minimized using MINUIT [18] with two free parameters 
Ub and Ug where ris = e^x Nq§ x E{B^ i^u) with ee being the efficiency in 
the fit region, and N^^ the total number of BB events analysed. We assume 
the number of the charged and neutral BB pairs to be equal. 

Figure 4 shows the Mbc distributions of events in the AE signal region together 
with the fit results. We observe 12 (15) events for the muon (electron) mode in 
the signal region. The signal yield extracted from the fit is 4.1 ±3.1 events for 
the muon mode and — 1.8±3.3 events for the electron mode in the signal region. 
For the SM branching fractions, we expect 2.8±0.2 and (7.3±1.4) x 10~^ events 
for the muon mode and the electron mode, respectively. The significance of 
the signal in the muon mode is 1.3, which is defined as y21n(Zi^^^^^yZo) where 
>Cmax is the likelihood value for the best-fit signal yield and Cq is the likelihood 
value for no signal event. No excess of events is observed in the electron mode. 
Selection efficiencies, expected numbers of events for signal and background 
and fit results are summarized in Table 1. 
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Fig. 4. distributions for the events in the AE signal region, together with the 
fit results (dotted lines). The solid curves are the background contributions. The 
dashed curves are the signal contributions. The signal contribution in the electron 
mode is multiplied by a factor of —4 to make it visible on the plot. 

5 Systematic Uncertainties 



The main sources of systematic uncertainty in calculating B{B^ l^v) arise 
from the uncertainties in the number of B^B~ events (A'^^^ = (276.6 ±3.1) x 
10^), determination of the signal efficiency, and parameterization of the Mbc 
distribution for the signal and background. The uncertainty due to the number 
of BB pairs is 1.1%. The uncertainties from the signal efficiencies are: 1% due 
to the uncertainties in the track-finding efficiency for the signal, 4.4% due to 
the uncertainty in the muon ID efficiency and 1.1% due to the electron ID 
efficiency, 2.3 % (2.1 %) for the muon (electron) mode from the MC statistics. 
We calculate an efficiency correction factor, to account for differences between 
MC and data, by analyzing a control data sample of fully reconstructed 5"^ 
r)(*)^Ti+ decays, where we treat the pion as a signal lepton and the D^*^^ as 
the accompanying neutrino. The event topology is similar to that of the signal 
events as both are two-body decays. The companion B is reconstructed in 
the control data samples as in the signal sample. We compare the efficiencies 
of the control data sample and a corresponding MC sample and determine 
correction factors for the signal efficiencies. The correction factors obtained 
are 1.01 ± 0.04 for the muon mode and 1.13 ± 0.04 for the electron mode. 
We apply a correction only to the electron mode, while uncertainties on the 
correction factors contribute to the systematic uncertainties of both modes at 
the level of 3.6 %. 

The uncertainties related to the signal Mbc shapes are estimated by repeating 
the fits while varying the Crystal Ball function parameters by their uncertain- 
ties; their contribution is 6.5 % for the muon mode and 3.2 % for the electron 
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Sources 




Muon Mode 


Electron Mode 






1.1% 


1.1% 


Signal Efficiency 


Lepton ID 


4.4% 


1.1% 




Tracking 


1.0% 


1.0% 




MU Statistics 


2.. 6/0 


O 1 0/ 

z.i/c 






3.6% 


3.6% 


Afbc Shape 


Signal 


6.5% 


3.2% 




Background 


8.1% 


15.7% 


Total 




12.2% 


16.7% 



Table 2 

Summary of systematic uncertainties 



mode. Uncertainties due to the background Mbc shapes are estimated in a 
similar manner by varying the ARGUS function parameters and yield 8.1% 
for the muon mode and 15.7% for the electron mode. Table 2 summarizes the 
contributions to the systematic uncertainties. 



6 Limits on Branching Fractions 



Figure 5 shows the dependence of the likelihood functions on the branching 
fractions, C{B), before and after inclusion of the systematic uncertainties. To 
account for systematic uncertainties in the calculation of the upper limits for 
each mode, we convolve the likelihood function C{B) with a Gaussian distribu- 
tion, where the sigma of the Gaussian corresponds to the size of the systematic 
uncertainty in this mode. The 90 % confidence level for the upper limit on the 
branching fraction, iSgo, is defined by 0.9 = J^'° C{B)dB / C{B)dB. 

We obtain the following upper limits on the branching fractions at the 90 % 
confidence level: 



^ /i+z/^) < 1.7 X 10"^ (90%C.L.) (5) 
^ e+i/e) < 9.8 X 10"^ (90%C.L.), (6) 

including the effect of the systematic uncertainties. The expected sensitivities 
on the upper limits with the present dataset, computed using toy MC studies 
with a null signal hypothesis, are 1.0 x 10~^ for the muon mode and 1.1 x 10~^ 
for the electron mode. 
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Fig. 5. Likelihood function dependence on the branching fractions. The solid and 
dotted curves represent the likelihood functions without and with inclusion of sys- 
tematic uncertainties, respectively. The arrows indicate the upper limits on the 
branching fractions at 90% confidence level. 

7 Conclusion 



We have searched for the purely leptonic decays l^^^ii and e^Ue 

using data collected by the Belle detector at the KEKB e'^e" asymmetric- 
energy collider. We have found no evidence of signal in either decay mode. 
We set upper limits on the branching fractions: B{B'^ — > /U."'"z/^) < 1.7 x 10~^ 
and B{B^ — >■ e^Vg.) < 9.8 x 10~^ at 90 % confidence level. These limits are 
the most stringent to date and improve the previous published limits [6,7] by 
a factor of 4 for B^ f^'^'^fj. and 15 for B^ e^Ve- 
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